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Abstract 

We solve the following problem. Given a continuous complex-valued potential qi defined on a 
segment [—a, a] and let 92 be the potential of a Darboux transformed Schrodinger operator, that 
is q2 — — gi-f 2(y)^ where / is a nonvanishing solution of the equation Aif = — qi{x)) f = 0. 
Suppose a transmutation operator Ti is known such that AiTiu = Ti for any u £ [—a, a]. 

^ , Find an analogous transmutation operator for A2 = — 92(2;). 

' It is well known that the transmutation operators can be realized in the form of Volterra 

I integral operators with continuously differentiable kernels. Given a kernel Ki of the transmuta- 

. tion operator Ti we find the kernel K2 of T2 in a closed form in terms of Ki. As a corollary 

' interesting commutation relations between Ti and T2 are obtained which then are used in or- 

der to construct the transmutation operator for the one-dimensional Dirac system with a scalar 
potential. 

1 Introduction 



, Transmutation operators are a widely used tool in the theory of linear differential equations (see, e.g., 

^ ' [2) ID; [23] > Hi], [3S] and the recent review [53]). It is well known that under certain quite general 

conditions the transmutation operator transmuting the operator A = —-^ +q{x) into B = —-^ is a 
Volterra integral operator with good properties. Its kernel can be obtained as a solution of a certain 
Goursat problem for the Klein-Gordon equation with a variable coefficient. In the book [8] another 
approach to the transmutation was developed. It was shown that to every (regular) linear second- 
order ordinary differential operator L one can associate a linear space spanned on a so-called L-basis 
- an infinite family of functions {<Pfc}fc°^o such that Ltpk = for fc = 0, 1, L(pk = k(k — l)(pfc_2, for 
fc = 2, 3, . . . and all i^k satisfy certain prescribed initial conditions. Then the operator of transmutation 
was introduced as an operation transforming functions from one such linear space corresponding to 
a certain operator L to functions from another linear space corresponding to another operator M, 
and the transformation consists in substituting the L-basis with the M-basis preserving the same 
coefficients in the expansion. 

In [3] the relation between the transmutation operator in the form of a Volterra integral operator 
and as a transmutation of the bases was clarified. In particular, it was shown that the corresponding 
Volterra integral transmutation operator transmutes the powers of the independent variable into 
the elements ^pk of an appropriate L-basis. In order to assure such an important mapping property 
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of the transmutation operator a parametrized family of transmutation operators for the Schrodinger 
operator A was introduced. This parametrized family resulted to be a necessary and natural tool 
for solving the problem of construction of the transmutation operator for a Darboux transformed 
Schrodinger operator. The solution to this problem is the main aim of the present paper. 

We give an explicit representation for the kernel of the transmutation operator corresponding to 
the Darboux transformed potential in terms of the transmutation kernel for its superpartner (Theorem 
1191) . Moreover, this result leads to interesting commutation relations between the two transmutation 
operators (Corollarv l25l) which in their turn allow us to obtain a transmutation operator for the one- 
dimensional Dirac system with a scalar potential as well as to prove (Theorem [251) the main property 
of the transmutation operator under less restrictive conditions than it has been proved until now. 
Namely, we show that the considered transmutation operators transmute A into B for any continuous 
(and not necessarily differentiable) potential q admitting a nonvanishing, in general, complex-valued 
solution (and this condition is always fulfilled, e.g., when q is real valued). We give several examples 
of explicitly constructed kernels of transmutation operators. It is worth mentioning that in the 
literature there are very few explicit examples and even in the case when g is a constant such kernel 
was presented recently in [3]. The results of the present paper allow us to enlarge considerably the 
list of available examples and give a relatively simple tool for constructing Darboux related sequences 
of the transmutation kernels. 

In Section [2] we introduce two main objects, the transmutation operators and the systems of 
recursive integrals or the L-bases. Besides recalling some recent results from [3] on the parametrized 
family of the transmutation operators we prove (Theorem ^ a criterion for a function to be a kernel 
of a transmutation operator from the parametrized family as well as a result on the explicit form 
of the inverse transmutation operator (Theorem 1101) . In Section [3] we present the construction of 
the transmutations for the Darboux transformed operators and obtain the commutation relations for 
them. We use the relations to weaken the conditions on the potential under which the transmutation 
operators preserve their main property (Theorem I26|). The section ends with several examples of 
explicitly constructed transmutation kernels. In Section U the results of the preceding section are 
used in order to obtain the transmutation operator for the one-dimensional Dirac system with a 
scalar potential. Finally, Section [5] contains some concluding remarks. 

2 Transmutation operators for Sturm-Liouville equations and 
systems of recursive integrals 

2.1 Transmutation operators 

According to the definition given by Levitan |23j, let i? be a linear topological space, A and B be 
linear operators: E ^ E. Let Ei and E2 be closed subspaces of E. 

Definition 1 A linear invertible operator T defined on the whole E and acting from Ei to E2 is called 
a transmutation operator for the pair of operators A and B if it fulfills the following two conditions. 

1. Both the operator T and its inverse T"^ are continuous in E; 

2. The following operator equality is valid 

AT = TB (1) 

or which is the same 

A = TBT^K 

Our main interest concerns the situation when A — B — and g is a continuous 

complex- valued function. Hence for our purposes it will be sufficient to consider the functional space 
E = [a, b] with the topology of uniform convergence. For simplicity we will assume that the interval 
is symmetric with respect to the origin, thus E = C^[— a, a]. 



2 



An operator of transmutation for such A and B can be realized in the form (see, e.g., [23] and 
[25] ) of a Volterra integral operator 

Tu{x)=u{x)+ [ K{x,t)u{t)dt (2) 



where K{x,t) = and H is the unique solution of the Goursat problem 



2 ' 2 

32 



gjg^ ^ + (3) 

1 r 

H{u,0)^- q{s)ds, H{0,v)=0. (4) 

If the potential q is continuously differentiable, the kernel K itself is the solution of the Goursat 
problem 

1 r 

K{x,x)^- q{s)ds, K{x,-x) = 0. (6) 

2 Jo 

If the potential q is n times continuously differentiable, the kernel K{x, is n + 1 times continuously 
differentiable with respect to both independent variables (see [2F|). 

An important property of this transmutation operator consists in the way how it maps solutions 
of the equation 

v" + iu^v = (7) 

into solutions of the equation 

u" - q{x)u + uj^u^O (8) 
where w is a complex number. Denote by eo(iw, a;) the solution of ([5]) satisfying the initial conditions 

eo(iw,0) = 1 and 60(10;, 0) — iuj. (9) 

The subindex "0" indicates that the initial conditions correspond to the point x = and the letter 
"e" reminds us that the initial values coincide with the initial values of the function e'"^. 
The transmutation operator ([2]) maps e'"^ into eo(iw,x), 

eo(zoJ,x)=T[e^n (10) 

(see ^ Theorem 1.2.1]). 

Following [25] we introduce the following notations 

Kc{x, t-h) = h + K{x, t) + K{x, -t) + h {K{x, C) - K{x, ~0}dC 

where h is a complex number, and 

Ks{x, t; 00) = K{x, t) - K{x, ~t). 

Theorem 2 ([25]) Solutions c{uj, x; h) and s(w, x; 00) of equation ([8]) satisfying the initial conditions 

c(a;,0;/i) = 1, 4(lj, 0; /i) = /i (11) 

s(w, 0; 00) = 0, 4(w, 0; cx)) = 1 (12) 
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can be represented in the form 

c{uj, x; h) = cos uix + j Kc{x,t; h) cos cot dt (13) 

and 



Denote by 



sin ux f , \ sin ut 

s{u},x;oo)^ h / Ks[x,t;oo) at. (14) 

u Jr. to 



Tcu{x) = u{x) + / Kc{x,t; h)u{t)dt 



Tsu{x) = u{x) + / Ks{x,t; oo)u{t)dt 



and 



the corresponding integral operators. As was pointed out in [3 they are not, in general, transmuta- 
tions on the whole space C^[—a, a]. 

2.2 A parametrized family of transmutation operators 

In [3] we introduced the parametrized family of operators T^, G C, given by the integral expression 

Thu{x) = u{x) + f K{x,t;h)u{t)dt (15) 

J —X 

where 



K{x, t;h) = f^+ K{x, t) + ^ [ {K{x, s) - K (x, -s)) ds. 
2 2 



(16) 



They are related to operators Ts and Tc (with the parameter h in the kernel of the latter operator) 

by 

Th^T,P,+TsPo, (17) 

where Pef{x) = (/(x) + /(—a;)) /2 and Pof{x) — [f{x) — /(—a;)) /2 are projectors to even and 
odd functions, respectively. In this subsection we show that the operators T/j are transmutations, 
summarize their properties and later, in Theorem 1151 we show how they act on powers of x. 
Let us notice that K(a;,t;0) = K(x,t) and that the expression 

h /■* 

K(x, t; h) - K(x, -t] h) = K{x, t) - K{x, -t) + - / {K{x, s) - K{x, -s)) ds = K{x, t) - K{x, -t) 

2 J -t 

does not depend on h. Thus, it is possible to compute K(a;, i; h) for any /i by a given K(x, t] hi) for 
some particular value hi. We formulate this result as the following statement. 

Theorem 3 ([3]) The integral kernels K.{x,t;h) andii.{x,t;hi) are related by the expression 

K{x,t;h) ^'^--^ +K{x,t;hi)+'^^-!^ f {Kix,s;hi) ~K{x,-s;hi)) ds. (18) 
2 ^ Jt 

Using pi?)) let us check how the operators Th act on solutions of ([7]). 

Proposition 4 The operator T/j maps a solution v of an equation v" +Ll>'^v — 0, where uj is a complex 
number, into the solution u of the equation u" — q(x)u + uj'^u — with the following correspondence 
of initial values 

w(0)=w(0), u'{Q)=v'{Q) + hv{Q). 
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Proof. It follows from ([T71) and Tlieoreni[2]that the operator Th indeed maps a solution of v" +uj'^v = 
into a solution of u" — q{x)u + uPu = 0. It is clear from the definition ([T5|) that 

Thv{{)) = v{{)). (19) 

As for the derivative, we have 

(Thv)\x) ^v'{x)+ I K'^{x,t;h)v{t)dt + K{x,x;h)v{x) +K{x,-x;h)v{-x), 



therefore 

{Thv)'{0) = v'{0) + K{0, 0; h)u{0) + K{x, -x; h)v{0) ^ v'{0) + hv{0). (20) 



Remark 5 As can be seen from the proof, formulas (|19|) and (|20|) are valid for any function v € 
C^[-a,a]. 

We know that the kernel of the transmutation operator T is related to the solution of the Goursat 
problem ©-(HI)- We show that a similar result holds for the operators T/j. 

Theorem 6 In order for the function K(x, t; h) to be the kernel of a transmutation operator acting 
as described in Proposition^ it is necessary and sufficient that H{u,v; h) := K{u + v,u — v;h) be a 
solution of the Goursat problem 

qIqJ ' =qiu + v)Hiu,v;h), (21) 

H{u,0;h) = ^ + ^£q{s)ds, H{0,v;h) = ^. (22) 

// the potential q is continuously differentiable, the function K{x, t; h) itself should be the solution of 
the Goursat problem 



— niT.') K(t. t- h\ = — 



yQ^2 q{^)jK{x,t-h)^—K{x,t-h). (23) 

K{x,x;h)='^ + ^J^ q{s)ds, K{x,-x;h) = ^. (24) 

Proof. If some transmutation operator Th acts as described in Proposition |4l then on a common 
domain of definition it coincides with the constructed operator Th, because both operators are contin- 
uous and act identically on the dense subset of functions, namely, on linear combinations of solutions 
of equations ([7]). Hence it is sufficient to verify that the kernel K(x,i;/i) satisfies conditions of the 
theorem. 

Suppose q S C^[—a,a]. Then the kernel K{x,t) is twice continuously differentiable and we can 
check pS)) directly. In the following calculation K'^{x, ±x) means K'^{x, t)\^__^^ and K[{x, ±x) means 
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{dl-d^,)nx,t-,h)^{dl~d^,) 



{K{x,s) - K{x,~s)) ds 



j{x)K{x,t) + (^j^ {K(x,s) ~ K'^{x,-s)) ds + K{x,x) - K{x,-x)j + 
+ ^dt {K{x, t) - Kix, -t)) = qix)K{x, t) + ^q{x) + 



+ K - KA^, -s)) ds + K'^ix, x) - K'^ix, -x) + K',{x, t) + K',{x, -t) ^ 

= q{x)K{x, t) + -q{x) + 2 f {K'Jsi^, s) - K'J^{x, -s) + q{x){K{x, s) - K{x, -s))) ds+ 
+ K'^ {x, x) - K'^ (x, ~x) + K[{x, t) + K[{x, -t)J = q{x) (k{x, t) + ^+ {K{x, s) ~K{x,-s)) ds^ + 

+ 1 {Kix, X) - Kix, t) + Ki^. -) - Ki^. -t) + Ki^, .) - Ki^, + Ki^, t) + Ki^. ^t)) = 

= q{x) [k{x, + ^ + ^ {K{x, s) - K{x, -s)) ds^ = qix)K{x, t; h), 

where we have used jS]), dH) and the equahties K'^{x,x) + Kf{x,x) — dx{K{x,x)) — l/2q{x), 
K'^{x,—x) — K[{x,—x) = dx{K{x,—x)) — 0. Equahties (1^ are vahd by the definition (ITBl) of 
'K{x,t;h). Clearly, the function H(u,w;/i) = K(m ^ v,u — v;h) is also the solution of the Goursat 
problem ((2T|) - ((22)) . 

Suppose now that q £ C[—a, a]. Consider a sequence of functions {qn}neti, where (/„ G C^[—a, a], 
n G N and qn ^ q, n ^ oo uniformly on [—a, a]. Denote by Hn{u,v) the solutions of Goursat 
problems ©-(U)) with the potentials qn, and by H{u, v) the solution with the potential q. It is known 
(see, e.g., [IS]), that Hn{u,v) — > H{u,v), n -> oo uniformly. The same is valid for the functions 
Kn{x,t) = Hn{^Y^, ^Y^) and K{x,t). Denote by K„(x,t;/i) the integral kernels constructed from 
Kn{x, t) by means of (ITS)) and let H„(u, v; h) = K„(u + v,u — v;h). Then H„(m, v; h) — H(w, v; h) — 
K(7i + v,u — v; h), n — > cx) uniformly. Similarly to |25[ 136] . the function H„('u, v; h) is the solution of 
the Goursat problem (I^T]) - (P^ with the potential (/„ if and only if H„(u,z;;/i) satisfies the integral 
equation 

(7„(a + /3)H„(a,/3;/i)d/3. (25) 



H„(u,u;ft,) = - + - / qniy)dy+ da 

^ ^ Jo Jo Jo 



From the first part of the proof we know that the functions H„(u,w; h) satisfy (pS)) . Passing to the 
limit when n — >■ cx) in (|25]) . we obtain that ll{u,v;h) satisfies (|25|) with the potential (7 and hence 
satisfies the Goursat problem (| 2 II) - ([22)) . 

For the second part of the theorem, suppose that a function H{u, v; h) satisfies the Goursat 
problem pi p -(|22 )) . Consider the Goursat problem (|3])-(|4]). It is well posed, and from its solution 
using (fTB| we can construct the kernel K.{x,t;h). From the first part of the proof we have that 
the function H{u,v,h) = K(u + v,u — v; h) satisfies the Goursat problem (|^T)) - (|2^ . and from the 
uniqueness of the solution it follows that H = H and H gives us the transmutation operator acting 
as described in Proposition U) ■ 

Example 7 Consider a function k{x,t) — 2{x+i) (^'^^^^j *^ Subsection \3.3\ it is explained how it can 
be obtained). We have 

t-1 2 t-1 



{dt-dt)k{x,t) 



{x + lf (:e + 1)2 ■ 2(:e + 1) ' 



-x-i 1 ^ , x-i I I r 2 

^("'^"^^ 2(^ = ^2 '^"'^ '=^"'")^2(^ = -2 + 2X (7TTF ' 
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thus the function k{x,t) satisfies the Goursat problem (j23p - (|24p with q(x) = 2/{x + 1)^ and h = — 1 
and by Theorem\^is the kernel of the transmutation operator T_i. 

Consider two solutions of the equation u" = 0, wi = 1 and U2 ~ x. We have that their images 

■yi = T_iui = 1 + / — -dt 



_^2{x + l) x + 1 
r (t-l)t , x^+3x^ + 5x 

V2 = T_iM2 =X+ '— dt - 



_^2{x+l) 3(a; + l) ' 

are solutions of the equation v" — 2/{x + l)'^v = with the initial values in agreement with Proposition 
For the solution U3 = cos a; of the equation u" = —u we get that 

rx^ r it— I) cost , sinx 
W3 = T_i-U3 = cosa; + / — — — —dt = cos x 



2{x+l) x + V 

is a solution of the equation v" — j^^.'^iyi v — —v. 

For the case q ^ [—a, a] the Voherra-type integral operator is a transmutation in the sense 
of Definition [1] on the space C^[— a, a] if and only if the integral kernel K{x,t) satisfies the Goursat 
problem 

K(x,t) ^ —K{x,t), 

1 

K{x,x)^C+- qis)ds, K{x,-x)=C, 
2 Jo 

where C is a constant (see [24l [35]). The proof of this fact consists in differentiating of the integral 
and integrating by parts, similar to the proof of Theorem [6l Due to Theorem [6] we see that T/i is 
indeed a transmutation. 

Theorem 8 Let q G C^[— a, a]. Then the operator T/i given by (jlSp satisfies the equality 



2 



^+q{x)]Th[u] = T 



dx 



■2 



d' 



(26) 



for any u £ C [— a,a 



Remark 9 This theorem was proved in under the additional assumption that there exists a non- 
vanishing solution f of f" — qf = on [—a, a]. Later, in Theorem\2^ we show that such additional 
assumption allows one to weaken the requirement on the potential q to q <E C[—a,a] keeping valid 
equality i26\) . 



As can be seen from (|15p . to define the transmutation operator we have to know its integral kernel 
in the region —a < a; < a, |t| < The corresponding function H should be defined in two domains, 
upper-right and lower-left triangles Hi: u > 0,v > 0,u -\- v < a and TI3 : u < 0,v < 0,\u\ -\- \v\ < a 
meanwhile the Goursat problem (|21|) -(|22 |) may be considered in four domains, triangles Hi and II3 
and triangles 112 : u < 0, w > 0, |w| -I- |u| < a, 114 : u > 0,v < 0, |m| -I- \v\ < a. The potential 
q is continuous in each of these domains, boundary conditions (|22p are continuously differentiable, 
therefore the corresponding solutions of the Goursat problems (see [3S]) belong to (7^(11^), i — 1 . . .A. 
Hence the function H(it, w;/i) is defined in the combined domain |m| -|- I?;! < a and belongs to the 
class on it. Consequently, the function K(a;,t;/i) — H(^i^,^^;/i) is defined on the domain 
n : —a < X < a, —a < t < a and is continuously differentiable there. For the rest of this paper we 
assume that the kernel K is defined on this larger domain 11. 
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Let us find the inverse operator T^^ ^ . Since is the Volterra integral operator its inverse (see 
[24| ) is again a Volterra integral operator 

Tl^u{x) ^u{x)+ I L{x,t;h)u(t)dt, 



where L{x,t;h) satisfies (for q e C^[~a,a]) the Goursat problem 

52 _ 



1 r 

L{x,x;h) = Ci — — I q{s)ds, L{x, 

2 Jo 



-x;/i) = Ci. (28) 



To determine the value of the constant Ci we use Proposition 01 Namely, we observe that the 
integral operator with the kernel L satisfying the Goursat problem (P7)) ~ (P5)) gives us the following 
correspondence of initial values 

w(0)=m(0), u'(0) = m'(0) +2Ciu(0), 

where v — T^^w. Therefore for the operator to be the inverse of the operator Th we should take 
Ci = -h/2. Comparing the Goursat problem (gTll-dlHl) with Ci = -h/2 with (gSll-dlll) we conclude 
that L{x, t; h) := — K(t, x; h) is the unique solution of (|27P " (p8)) and hence it is the kernel of the inverse 
operator. The assumption q G [—a, a] is not essential and can be easily overcome by considering an 
approximating sequence of continuously differentiable potentials g„. Thus, the following statement is 
proved. 

Theorem 10 The inverse operator can be represented as the Volterra integral operator 

Tl^u{x)^u{x)~ [ K{t,x;h)u{t)dt. (29) 



Example 11 For the operator T_i with the kernel K{x,t) and functions Ui and Vi, i = 1,2, 3 from 
Example it is easy to see that indeed 

u,{x)^TZ\[v,]{x)=v^(x)- j -^Z^^y^[t)dt, i = l,2,3. 
2.3 A complete system of recursive integrals 

Let / S C^{a, b) D C^[a, b] be a complex valued function and f{x) 7^ for any x E [a, b]. The interval 
(a, 6) is supposed to be finite. Let us consider the following auxiliary functions 

X(o)(a;) = X(o)(a;) = 1, (30) 

X^^\x)=n rx("'i)(s)(/2(s))^-'^"" ds, (31) 

J Xo 

X(")(x) =n rx("-i)(s)(f(s))'"''" ds, (32) 



where Xq is an arbitrary fixed point in [a,b]. We introduce the infinite system of functions {^k}'kLo 
defined as follows 

ff{x)X(''\x), fcodd, 
\f{x)X^''>{x), k even, 

where the definition of X^''^ and X'*^^ is given by (p0l) -(|32 |) with xq being an arbitrary point of the 
interval [a,b]. 
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Example 12 Let f — 1, a = 0, b = 1. Then it is easy to see that choosing xq — we have 
fk{x) — , k G No where by Nq we denote the set of non-negative integers. 

In |19j it was shown that the system {'Pk}'kLo is complete in ^2(0, b) and in |20| its completeness 
in the space of piecewise differentiable functions with respect to the maximum norm was obtained 
and series expansions in terms of the functions ipk were studied. 

The system (|33p is closely related to the notion of the L-basis introduced and studied in [8 . Here 
the letter L corresponds to a linear ordinary differential operator. This becomes more transparent 
from the following result obtained in [T7] (for additional details and simpler proof see [TB] and [?!] ) 
establishing the relation of the system of functions {^k}k=a to Sturm-Liouville equations. 

Theorem 13 ([17J) Let q be a continuous complex valued function of an independent real variable 
a; e [a,b], A be an arbitrary complex number. Suppose there exists a solution f of the equation 

f" -qf = (34) 

on (a^b) such that f G C^(a,b) together with 1/f are bounded on [a,b]. Then the general solution of 
the equation 

u" — qu = \u (35) 

on (a, b) has the form 

U = CiUi + 02^2 

where ci and C2 are arbitrary complex constants, 

"i = E(^^2. and ^2 = E(^^^^2fe+i (36) 



k=0 ' ' k=0 

and both series converge uniformly on [a, b] . 

Remark 14 It is easy to see that by definition the solutions ui and U2 satisfy the following initial 
conditions 

ui{xa) = /(xo), u\{x()) = /'(a;o), (37) 
M2(a;o) = 0, U2(a;o) = l//(a;o). (38) 

2.4 Transmutations and systems of recursive integrals 

In this subsection we show the connection between the transmutation operators T/j and the functions 
Lpk. All the results of this subsection were proved in [3]. Here we recall them with explanations. 

We suppose that / is a solution of ([M)) fulfilling the condition of Theorem [T3] on a finite interval 
(—a, a). We normalize / in such a way that /(O) — 1, and let /'(O) = h where h is some complex 
constant. Let us obtain the expansion of the solution c(w,a;; h) from Subsection 12.11 in terms of the 
functions ipk- According to Remark [14] the solutions ((36|) of equation (|35|) have the following initial 
values 

Mi(0) = 1, m;(0) = ft, M2(0) = 0, 7/2(0) = 1. 
Hence due to (fTTj) we obtain c(aj, x; h) — ui{x). From (fT3|) and ([36]) we have the equality 

^ ^-^H^2kix) = V ^ + / K,(x, t; h) V ^ dt. 



As the series under the sign of integral converges uniformly and the kernel Kc{x,t;h) is at least 
continuously differentiable (for a continuous q |25j) we obtain the following relation 
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This equality holds for any u hence we obtain the termwise relations 

V2k=T,[x''% fee No. (39) 

Similarly, due to (|12p we observe that s(w,x;cxd) = 1*2(0;) and from (|14p we find that 

(P2fc+i =T,[x2^+i], fee No. (40) 

Hence from the last two equalities and pT|) we conclude that the following statement is true. 

Theorem 15 ([3j) Let q he a continuous complex valued function of an independent real variable 
X £ [—a, a], and f he a particular solution of \34^ such that f G C^(— a, a) together with 1/f are 
bounded on [—a, a] and normalized as /(O) = 1, /'(O) — h, where h is a complex number. Then the 
operator U5\) with the kernel defined by U6]) transforms x^ into (pk{x) for any fc €E No. 

Thus, we clarified what is the result of application of the transmutation T/j to the powers of the 
independent variable. This is very useful due to the fact that as a rule the construction of the kernel 
K(2:,t;/i) in a more or less explicit form up to now is impossible. Our result offers an algorithm 
for transmuting functions which can be represented or at least approximated by finite or infinite 
polynomials in the situation when K(x,t; h) is unknown. 

Remark 16 Let f be the solution of ^34^ satisfying the initial conditions 

/(0) = 1 and /'(0)=0. (41) 

// it does not vanish on [—a, a] then from Theorem 1 1 5\ we obtain that the original transmutation 
operator T transmutes powers of the independent variable into (pk- fk — Tlx'^] for any fc G No. 
In general, of course f may have zeros on [—a, a] and hence one can guarantee that the operator 
T transmutes the powers of x into ipk whose construction is based on the solution f satisfying |^_?[ ) 
only in some neighborhood of the origin. If we construct (pk by a solution f such that /(O) 1 
and /'(O) = h, then the original transmutation operator T does not map powers of the independent 
variable exactly to the functions (fk. As it was shown in fS^ the following equalities are valid 

Lfk ~ T[x^] when k is odd, 

fk — -, Tfk+i = Tlx'^] when fc G Nq is even. 

fc + 1 

This observation explains the necessity to consider the parametrized family of transmutation operators 
whose mapping properties are well adjusted to corresponding families of the functions Lp^ ■ 

Example 17 For the operator T_i from Example^ consider the function f = T_i[l] = as a 
solution of (p4|) such that /(O) = 1 and /'(O) = h — —1, nonvanishing on any [—a, a] C (—1, 1). The 
first 3 functions (fk are given by 

1 x^ + _^ 3^ 2x^ + 3a;2 

'^0 = /-^' 3(x + l) ' 3(x + l) ■ 

We have seen already that fQ — T_i[l] and fi ~ T_i[a;]. We have also T_i[a;2] — x'^ — ^(^^^y-^ = f2- 

3 Transmutation operators for Darboux transformed equa- 
tions 

3.1 The Darboux transformation 

,2 

Consider a Sturm-Liouville operator Ai := ^ where qi is a continuous complex-valued 

function on some segment [—a, a]. Suppose a solution / of the equation Aif = is given such that 
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/(.t) 7^ 0, X G [—a, a], it is normalized as /(O) = 1 and h /'(O) is some complex number. Then the 

-qi{x), is known as the Darboux transformation 



■q2ix), where q2ix) = 2(^^; 



operator A2 := - 
of the operator Ai. 

Initially the Darboux transformation served for establishing a relation between the general solution 
u of the equation Aiu = Xu and the general solution v of the equation A2V — Xv hy a known particular 
solution / of the equation Aiu = 0. This relation is given by the formula v{x) — u'{x) — u{x)y^. 
Later on it was found that the Darboux transformation was closely related to the factorization of the 
Schrodinger equation, and nowadays it is used in dozens of works, see e.g. [71 HUl ESI |3^ in connection 
with solitons and integrable systems, e.g. [TJ [11] [28l [29] and the review [3T] of applications to quantum 
mechanics. 

We remind some well known facts about the Darboux transformation. First, 1// is the solution 
of A2U = 0. Second, it is closely related to the factorization of Sturm-Liouville and one-dimensional 
Schrodinger operators. Namely, we have 



A, 



f 



/ ^ / 



(42) 
(43) 



Suppose that u is a solution of the equation Aiu = lou for some w G C. Then the function v 
ipx — ^)'^ — {f9xj)u is a solution of the equation A2V — ujv, and vice versa, given a solution v 
A2V — (jjv, the function u = [dx + y)w = {jdxf)v is a solution of Aiu = uju. 



3.2 Construction of transmutation operators for Darboux transformed 
equations 

Consider the operators Ai, A2 and a particular solution / of the equation Aif — as in the previous 
subsection. Suppose that the operator Ti-h which transmutes the operator Ai into the operator 
B = (P/dx^ is known in the sense that its kernel 'Ki{x,t; h) is given. As before h — /'(O) and T^h 
transforms solutions according to Proposition [Jj As the parameter h is fixed by the value /'(O) we 
will simply write Ti instead of Ti^^,. We know that the function 1// is the non- vanishing solution of 
the equation A2U — satisfying l//(0) = 1 and (l//)'(0) = —h. Hence it is natural to look for the 
operator T2:~h transmuting the operator A2 into the operator B. We will simply write T2 further 
in this paper. 

Let us explain the idea for obtaining the operator T2 . We want to find an operator transforming 
solutions of the equation Bu + lu'^u = into solutions of the equation A2U + ui'^u = 0, see the first 
diagram below. Starting with a solution <t of the equation (9^ +i^^)(T = 0, by application of Ti we get 
a solution of {Ai +uj'^)u — 0, and the expression {fdxj)Ti<J is a solution of {A2 +uj'^)v ~ 0. But the 
operator (^fdxj)Ti is unbounded and hence cannot coincide with the operator T2. In order to find 
the required bounded operator we may consider the second copy of the equation (9^+a;^)u = 0, which 
is a result of the Darboux transformation applied to {d^ + ijj^)a = with respect to the particular 
solution 5=1 and construct the operator T2 by making the second diagram commutative. In order 
to obtain a bounded operator T2, instead of using fdxj for the last step, we will use the inverse of 
i.e. i(/;/(s). d5 + C). 



Ti 



91 



dl-q2 



dl-qi 



12 



I , .2 
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That explains how to obtain the following theorem. 
Theorem 18 The operator T2, acting on solutions u of equations {d^ +lu'^)u — 0, cu CI C by the rule 

T2[u]{x) f[rj)T,[u'Md7j + uiO)^ (44) 

coincides with the transmutation operator T2;-h- 

Proof. According to Proposition |4] the operator T2:-/i transforms a solution u of the equation 
(5^ + uj'^)u = into a solution of the equation {A2 + uP')v = {d^ — q2{x) + lu'^)v = with the 
correspondence of initial conditions 

T2;-/,M(0) = u(0), (T2;-ftM)'(0) = u'iO) - hu{0). (45) 

Hence we have to verify the same properties for the operator T2 . Consider a solution u of the equation 
(9^ + ui^)u = for some particular w € C. The function u' is again a solution of this equation, and 
by Proposition m Tiu is a solution of the equation {d^ — Qi + uP')v = 0, so with the aid of the 
factorization ([^^ we obtain 

j9./'9.jTiM = -c.2TiK]. (46) 



Let us apply the operator A2 to the function T2U and use the factorization (|43p . 

A2T2[u]{x) = [fd^^dj) f{s)T,[u']{s)ds + u{Q)yj = (/5,i)(TiM). (47) 

Applying jdxj to both sides of (ITTl) due to PS|) we obtain 

Hence the function A2T2[u] may differ from the function —u!'^T2[u] only by c//, where c is a constant. 
To find the value of the constant we compute the values of both expressions for a; = 0. We have 



T2M(0) =u(0) 

and 

A2T2[u]{x) = (/a i)(TiM(x)) = -IMt,[u']{x) + (TiM(x))' 



fJ^ - - fix) 

As u' is a solution of (d^ + io'^)v = 0, by Proposition H we have Ti[m'](0) = u'(0), (Ti[u'])'(0) = 
m"(0) + hu'{0) = —uj'^u{0) + hu'{0) (the last equality holds since the function u is also a solution of 
{d^+uj^)v = 0). Therefore 

^272 M(0) = -W(0) + wii(0) + W(0) = -^2^(0). 

Hence c = and the operator T2 maps solutions of the equation {d^ + uj^)u = into solutions of the 
equation {A2 +uj'^)v = for any u! G C To finish the proof, we have to check conditions (|45|) for the 
operator T2. We have 

T2M(0)^^=u(0) 

and 

(T2M)'(x) = -^T2M(x) + • /(x)TiM(x) = -/Mt2M(x) + T,[u']{x), 
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hence for x — we obtain 

(T2M)'(0) = ~hT2[u]{0) + Ti[u']{0) = ~hu{0) + u'{0), 

which finishes the proof. ■ 

Now we show that the operator T2 can be written as a Volterra integral operator and, as a 
consequence, extended by continuity to a wider class of functions. 

Theorem 19 The operator T2 admits a representation as the Volterra integral operator 

T2[u]{x)^u(x)+ [ K2{x,t;-h)u{t)dt, (48) 



with the kernel 

K2(x, t;-h) = ^j^(^ J^^ dtK^is, t- h)f{s) ds + ^/(-i)) ■ (49) 
Proof. Consider the expression 

f{t)Ti[u']{t)dt + u{0) = f{t)u'{t)dt+ f{t) Ki{t,s;h)u'{s)dsdt + u{0). 



Suppose that x > (the opposite case is similar) . We integrate by parts the first integral and change 
the order of integration in the second integral, 

f{t)Ti[u']{t)dt + u{0) = 

nX nX pX 

f{x)u{x)~ f{0)u{0)- f'{t)u{t)dt+ u'{s) Ki{t,s;h)f{t)dtds + u{0) = 



f{x)u{x)- / f'{t)u{t)dt + u{s) / Ki{t,s;h)f{t)dt 



- f I Ki (t, s; h)f{t) dt] ds = 

J-x ds\J^^\ J 

f{x)u{x)-J^ f{t)u{t)dt-J uis)^(^J Ki{t,s;h)f{t)dtys. (50) 



To continue the proof we will use the extension of the kernel Ki onto the square \x\ < a, \t\ < a 
as a continuously differentiable function. Consider the integral 

r\s\ 

j Ki{t,s;h)f{t)dt. 

For s < it equals zero. For s > note that by Theorem[Tni Ti[l](a:) = f{x), so we obtain from 

l = T^^[f]{x)^ fix) - r K^{t,x-h)fit)dt. 

Hence 

"I**' f if s < 

K^{t,s;h)f{t)dt^{ : - ' (51) 

/(s) - 1, if s > 0. 
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Combining ([50]) and we get 



K,it,s;h)f{t)dt 



/(t)TiK](t) + u(0) = f{x)u{x) - I f'{t)u{t) dt- 

y' ^ Ki{t,s;h)f{t)d?jds 
= f{x)u{x)^l^ f'{t)u{t)dt- j <s)^^fyj Ki{t,s;h)f{t)dt^ds + u{s)f'{s)ds = 
^f{x)u{x)~j u{s)l^j ^K,it,s;h)f{t)dt + K,i-s,s;h)f{-s)ys = 

^f{x)u(x)^ j u{s)(^j ±K^{t,s-h)f{t)dt+^f{-s)^ds. 

Hence 



T2[u]{x)^j^\^j^ f{'n)Ti[u'MdTi + u{0)j =u{x)+ J K2{x,t;-h)u{t)dt, 
where the kernel K2 is given by (HH]) . ■ 

Remark 20 Note that in the proof of Theorem{I^ we made use only of the fact that u € C^[— a, a] 
and never required that u he a solution of the equation d'^u + uPu = 0. Therefore both representations 
for the operator T2 obtained in Theorems \18\ and \19\ coincide on any function u (1 C^[— a,a]. 



Remark 21 // the integral kernel Ki is known only in the domain \x\ < a, \t\ < \x\, from ()50p 
it is also possible to obtain the expression of the integral kernel for the Volterra integral operator 
representing T2 ■ Since 

^ Ki(s, t- h)f{s) ds^ - ^Ki(s, t- h)f{s) ds - Ki(|i|, i; h)f{\t\) ■ {\t\)' = 

^ { j^^i |Ki(s, t- h)f{s) ds - ^if{t) ~ ^ X)* q,{s) ds for t > 0, 
\ J^^i |Ki(s, t; h)f{s) ds + yi-t) for t < 0, 



the integral kernel for x > is given by the expression 
K2{x,t] -h) = 



-7M (^'(*) + ft (Ki (s, t; h))f{s) ds - - q,{s) ds) tft>0, 

-TRl/ltl i{KAs,t;h))f{s)ds+yi~t)) ^ft < 0. 



Corollary 22 The operator T2 given by ()48p with the kernel (|49p coincides with T2 on C[— a, a]. 

Proof. By Theorems [18] and [19] the Volterra operators T2 and T2 coincide on the set of finite linear 
combinations of solutions of the equations (9^ + oj^)?/ = 0, uj G C. Since this set is dense in C[—a, a] , 
by continuity of T2 and T2 we obtain that they coincide on the whole C[— a, a]. ■ 
The next corollary follows immediately from Remark 

Corollary 23 The operator T2 given by (|44l) coincides with T2 on C^[~a,a]. 

Operator Ai is the Darboux transformation of the operator A2 with respect to the solution 1//, 
hence we obtain another relation between the operators Ti and T2. 
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Corollary 24 For any function u G C^[-~a,a] the equality 

1 



Ti[u]{x) ^ fix) 



fiv) 



T2[u']irj) drj + u{0) 



(52) 



is valid. 



From the second commutative diagram at the beginning of this subsection we may deduce some 
commutation relations between the operators Ti, T2 and d/dx. The proof immediately follows from 
(m and dSll). 



Corollary 25 The following operator equalities hold on C^[— a,a] 

a,/T2 = fTA 



f f 



(53) 
(54) 



Commutation relations (|53|) and (|54|) allow us to prove a more general version of [Sj Theorem 11] 
which is also Theorem [8] with different conditions on the potential q. 



Theorem 26 Under the conditions of Theorem ] 15\ the operator ()15l) with the kernel defined by (fT6 
satisfies 

d^ \ 
— + g(x) ) Th[u]^Th 



dx'^ 



for any u £ C^[— a,a]. 



Proof. Let T2 be the transmutation operator of the Darboux transformation of the operator A = 
— + q{x) with respect to the solution /. Let u e C^[—a, a], then u' € C^[—a, a] and we obtain by 
dSf , (IMl) and (USD 



dx"^ 



I d r n 1 d , d 1 f d'^ , , , 



In [20] the following notion of generalized derivatives was introduced. Consider a function g 
assuming that both / and g possess the derivatives of all orders up to the order n on the segment 
[—a, a]. Then in [—a, a] the following generalized derivatives are defined 

7fc(5)(x) = (/2(x))(-'^"'(7.-i(5))'(x) 



for k = 1,2, ... ,n. 

Let a function u be defined by the equality 



and assume that u G C"'[— a,a]. Note that below we do not necessarily require that the functions / 
and g be from C"[— a,a]. With the use of ([55]) and (|54l) we have 



71(3) = /' • (7T1")' = /' • jT^u' = fT2u', 

72(g) 



1 



f 

fT2u' 



1 



fT,u" = -Tiu". 



By induction we obtain the following corollary. 
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Corollary 27 Let u e C"[-a,a] and g = jTiu. Then 

j^(g) = /T2u'^''^ i/ k is odd, k <n, 

and ^ 

lk{g) — -jTiu^'^'^ if k is even, k < n. 

To finish this subsection, let us consider how the operator T2;^h acts on powers of the independent 
variable. By Theorem [T5l we have to construct the system of functions {fk} by formulas (l30 |) -([32 |) 
and ([33|) starting with the function 1//. As it can be seen from (p0|) ~([32 |) . for any nonnegative integer 
k we have 

X[^\^xf and X[f^.=xf\ 

where the subindex / or 1/f corresponds to the starting function used in ((30|) -p2 | . Therefore 
constructing the functions X^''\ X^^'^ according to ((30l) - ((32)) and defining 

r 1 XW, fc even, 

(here the "second half of the formal powers (|30 |) -([32 l) is used, cf. ([33])) we obtain the following 
statement. 

Proposition 28 Let a potential qi and a function f be as in Theorem 1 1 5i Let the operator T2--h 
be the transmutation operator for the Darboux transformed operator A2. Then T2;-h transforms 
into 4'k{x) for any fc G Nq. 

3.3 Examples 

We start with the operator — d^ /dx"^ . We have to pick up such a solution / of the equation 
^o/ = that /'// 7^ 0. This is in order to obtain an operator Ai 7^ as a result of the Darboux 
transformation of Aq. As such solution consider, e.g., /o(a;) = x + Both /g and I//0 are bounded on 
any segment [—a, a] C (—1; 1) and the Darboux transformed operator has the form Ai — — (x+i)'^ • 
The transmutation operator T for Aq is obviously an identity operator and Kn{x, t; 0) = 0. Since 
/q(0) = 1, we look for the parametrized operator Toa- Its kernel is given by (|18p : 'Ko{x,t; 1) = 1/2. 
From Theorem 1191 we obtain the transmutation kernel for the operator Ai 

K.(.,,;-l).--l^-l^ = J^. M 

the kernel from Example [T) 

To obtain a less trivial example consider again the operator Ai = — i^^^iy and the function 
fi{x) = {x + lY as a solution of Aif = 0. Since h = /{(O) = 2, we compute Ki(a;,i;2) from ([56)) 
using ([HI) 

3^2 + 6a; + 4 - it^ + 2t 
Ki(x,i;2 = . 

4(.T + 1) 

The Darboux transformation of the operator Ai with respect to the solution /i is the operator 
A2 = — j^x+i)^ ^^'^ Theorem [Tni the transmutation operator T2:-2 for A2 is given by the 
Volterra integral operator ([T5|) with the kernel 

K2(x f -2) = ^ f r ^^is + I? ds + (l- t)A = i^t-l){x + ir-3{t~mt + l) 
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This procedure may be continued iteratively. Consider operators 

n{n + 1) 



da;2 (x + 1)2 



/„(x) = {x + 1)"+^ is a solution of the equation Anf — 0. The Darboux transformation of the 
operator An with respect to the solution /„ is the operator 

d^_ (mx)_Y,ri{n + l)_d' {n+l){n + 2) 
dx^ \Ux)) 



2 



U{x)J (a; + l)2 dx^ (x + 1) 

i.e., exactly the operator A^+i- If we know K„(a;,t;— n) for the operator An, by we compute 
the kernel K„(x,i;n + 1) corresponding to the solution /„(x) and by Theorem \W\ we may calculate 
the kernel K„_|_i(x,t; —n — 1). Careful analysis shows that we have to integrate only polynomials in 
all integrals involved, so the described procedure can be performed up to any fixed n. 
Consider the Schrodinger equation 

u" + 2 sech^ {x)u = u. (57) 

This equation appears in soliton theory and as an example of a reflectionless potential in the one- 
dimensional quantum scattering theory (see, e.g. [22]). Equation (|57| can be obtained as a result of 
the Darboux transformation of the equation u" = u with respect to the solution f{x) = coshx. The 
transmutation operator for the operator Ai — &l — \ was calculated in [S', Example 3]. Its kernel is 
given by the expression 

l ./x^h{Vx^) 

= —2 ' 

where Ii is the modified Bessel function of the first kind. Hence from Theorem [19] we obtain the 
transmutation kernel for the operator yl2 = ^2 + 2 sech^ a; — 1 



^^^"'^'O) = 2^^) J_, [ s-t + ) 

4 Transmutation operator for the one-dimensional Dirac equa- 
tion with a Lorentz scalar potential 

One-dimensional Dirac equations with Lorentz scalar potentials are widely studied (see, for example, 
[5] [S] [H [H [H [13 [H [^ [33 [35 and for intertwining techniques for it) . 

According to |27| the Dirac equation in one space dimension with a Lorentz scalar potential can 
be written as 

{d,+m + Six))^i^Ei:2, (58) 
{-d^+m + S{x))yj2^E^Pi, (59) 

where m (m > 0) is the mass and S{x) is a Lorentz scalar. Denote rj = m + S and write the system 
(|59p in a matrix form as 



In order to apply the results on the transmutation operators and factorizations ([42]) . (|43| we consider 
a function / such that 

fix) , 
— — = -77 = -m- S{x). 

fix) 
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We can take f{x) — exp {m + S{s)) ds) , then /(O) = 1 and / does not vanish. Suppose the 



operators Ti and T2 are transmutations for the operators Ai — {d^ + y) {px — y) and A2 



{px — y) {px + y) respectively (corresponding to functions / and 1// in the sense of Proposition 
We will look for a solution of equation (IMl) in the form 

^ /Ti \ (ui 

where ui and U2 are some functions. From the commutation relations (j53p and (|54p we have 




1\ /Ti \ 

[ T2 U2 



/T"^ \ 

By multiplying both sides by the inverse matrix ^ t^^ y obtain 
dx 0\fui\ ^/T7^ OWO 1\ fTi 0\ fui 



E 



Therefore the operator ) transmutes any solution ( ) of the system 

u'l = £;u2 (61) 
u'2 = -Eui (62) 

into the solution ( ) of the system (15^1) with the initial conditions V'i(O) = ui(0), "02(0) = 1*2(0). 
And vice versa if ( ) is a solution of the system (l58t . (|59| . then the operator ^ t" f'"^'^'^™'^tes 
it into the solution ("^ ) of ^ such that lii(O) = 0i(O), ^2(0) = 02 (0). 



5 Conclusions 

An explicit representation for the transmutation operator corresponding to a Darboux transformed 
Schrodinger operator is given in terms of the transmutation kernel for its superpartner and as a 
corollary the transmutation operator for the one-dimensional Dirac system with a scalar potential is 
obtained. Several examples of explicitly constructed transmutation operators are given. We expect 
that the techniques developed in the present paper will be used in practical applications of the 
transmutation operators. 
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